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We study cold dense quark matter and hybrid neutron stars with a Dyson-Schwinger quark model and various 
choices of the quark-gluon vertex. We obtain the equation of state of quark matter in beta equilibrium and 
investigate the hadron-quark phase transition in combination with a hadronic equation of state derived within 
the Brueckner-Hartree-Fock many-body theory. Comparing with the results for quark matter with the rainbow 
approximation, the Ball-Chiu ansatz and the 1BC ansatz for the quark-gluon vertex lead to a reduction of the 
effective interaction at finite chemical potential, qualitatively similar to the effect of our gluon propagator. We 
find that the phase transition and the equation of state of the quark or mixed phase and consequently the resulting 
hybrid star mass and radius depend mainly on a global reduction of the effective interaction due to effects of 
both the quark-gluon vertex and gluon propagator, but are not sensitive to the vertex ansatz. 
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I. INTRODUCTION 

The possible appearance of quark matter (QM) in the inte¬ 
rior of massive neutron stars (NS) is one of the main issues 
in the physics of compact stars (U. Recent observations con¬ 
firm the existence of two NS of about two solar masses Hi. 
Based on a microscopic nucleonic equation of state (EOS), 
one expects that in such heavy NS the central particle density 
reaches values larger than 1 /fm 3 , where in fact quark degrees 
of freedom are expected to appear at a macroscopic level. 

The mass of a NS can be calculated by solving the Tolman- 
Oppenheimer-Volkoff (TOV) equations with the relevant EOS 
as input, which embodies the theoretical information of our 
theory on dense matter. The hybrid EOS including both 
hadronic matter and QM is usually obtained by combining 
EOSs of hadronic matter and QM within individual theo¬ 
ries/models. Unfortunately, while the microscopic theory of 
the nucleonic EOS has reached a high degree of sophistication 
Hill , the QM EOS is still poorly known at zero tempera¬ 
ture and at the high baryonic density appropriate for NS. 

Continuing a set of investigations using different quark 
models HHl, in Ref. JH we developed a Dyson-Schwinger 
model (DSM) for QM based on the Dyson-Schwinger equa¬ 
tions (DSE) of QCD JTMll. In that work we used the 
‘rainbow’ approximation, i.e., the bare quark-gluon vertex 
was employed. In combination with a baryonic EOS devel¬ 
oped within the Brueckner-Hartree-Fock (BHF) many-body 
approach of nuclear matter, we found that hybrid NS of two 
solar masses could be obtained when hyperons were not in¬ 
cluded and the nuclear matter EOS was stiff enough. 

Though the DSE are well based on QCD, in practice one 
has to work within certain truncation schemes of the infinite 
hierarchy of coupled equations. In our truncation scheme for 
the DSE of the quark propagator, we have thus to use an ansatz 
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for the gluon propagator and the quark-gluon vertex. In this 
paper, we investigate the effects of different choices for the 
quark-gluon vertex, employing the so-called Ball-Chiu (BC) 
vertex |2U l or the 1BC vertex, i.e., the first term of BC vertex. 
The comparison with the rainbow approximation will also be 
made. 

The paper is organized as follows. In section lHAl we briefly 
discuss the baryonic EOS in the BHF approach. Section HT bI 
concerns the DSM with the BC and 1BC ansatz. In section 
nm we present the results regarding NS structure, connecting 
the baryonic and QM EOS for beta-stable nuclear matter with 
a phase transition under the Gibbs construction. Section ||V] 
contains our conclusions. 


II. FORMALISM 

A. Hadronic matter within Brueckner theory 

Our EOS of hadronic matter obtained within the BHF ap¬ 
proach J2j]] has been amply discussed in previous publica¬ 
tions M- The basic input quantities of the calculation are the 
nucleon-nucleon two-body potentials, namely Argonne V\s 
ll22h . Bonn B l23ll . or Nijmegen 93 lj24j], supplemented with 
compatible three-body forces HUlIm 

This approach has also been extended with inclusion of hy¬ 
perons H [27| [281], which might appear in the core of a NS. 
The hyperonic EOS in this theory is very soft, which results 
in too low maximum masses of NS |29li] and often suppresses 
the appearance of quark matter. In this work we do not dis¬ 
cuss this aspect, but choose a purely nucleonic EOS with the 
Bonn B potential, supplemented by a compatible microscopic 
three-body force Hi- Since our main goal is to investigate the 
influence of our quark model, and the inclusion of quark mat¬ 
ter usually decreases the maximum mass of NS, we select a 
hard EOS of hadronic matter, which supports a large NS max¬ 
imum mass. 
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The BHF calculations provide the energy density e of the 
bulk system as a function of the relevant partial densities p„ 
from which all other thermodynamical quantities can be ob¬ 
tained, in particular chemical potentials and pressure. 



(1) 


2 d e de 


p(Pb)=Pb-7——=Pb-7—-£ = PbI^b-£- (2) 
dpB Pb dps 


The parameterized energy density function can be found in 
Ref.ll. 


B. Quark phase with the Dyson-Schwinger model 

For the deconfined quark phase, we adopt models based on 
the DSE of the quark propagator. At finite chemical potential 
p=p q = Pb/ 3, the quark propagator assumes a general form 
with rotational covariance, 

S{p;p)~ l = [nP + iMP4 + ip) + m q ] + I.(p-,p) (3) 

= ijpA(p 2 .p- u) + B(p 2 ,p ■ u) 

+iY 4 (P4 + ip)C(p 2 ,p-u), (4) 

where m q is the current quark mass, u = (0 ,ip), and possibil¬ 
ities of other structures, e.g., color superconductivity 1301431 . 
are neglected. The quark self-energy can be obtained from the 
gap equation, 

f d^ 

£(p;M) = J S 2 {p)D P c{p-q\p) 

X a ?i a 

x— YpS(q;p)—r a {q,p-,p), (5) 

where X a are the Gell-Mann matrices, g(p) is the cou¬ 
pling strength, D pa (k-,p) the dressed gluon propagator, and 
r a (q,p-,p) the dressed quark-gluon vertex at finite chemical 
potential. 

This vertex is a complicated quantity, which at zero chem¬ 
ical potential can be decomposed into 12 linearly indepen¬ 
dent Lorentz covariants. The DSE of the quark-gluon ver¬ 
tex depends on higher-order Green functions, and so far lit¬ 
tle is known about it. Commonly, one develops an ansatz 
constrained by all available reliable information, e.g., the 
Slavnov-Taylor identities (STI), effects of dynamical chiral 
symmetry breaking, and some recent lattice results I33l437tl . 
In our previous work 0 we employed the simple ‘rainbow’ 
approximation, T a (q,p-,p) = 

In this work, we will investigate another popular model, the 
so-called BC vertex ll20tl . The BC vertex satisfies the Ward- 
Takahashi identity (WTI) of QED, which constrains the lon¬ 
gitudinal part of the fermion-photon vertex, and is free of ki¬ 
netic singularity. Furthermore, the scalar part of the BC vertex 
gives a direct representation of the dynamical chiral symme¬ 
try breaking in the vertex. Though the WTI is different from 
the STI of QCD, this ansatz for the quark-gluon vertex is of¬ 
ten used as a starting point to constrain the vertex with STI of 


QCD [38]. Phenomenologically, it also improves the results 
of calculating hadron properties |[l9l.[39ll . It is widely used in 
comparison with the rainbow approximation to identify robust 
features of a vertex ansatz. 

The form of the BC vertex at finite chemical potential is 
developed in [4C|, 


ir l c (q,p',p) = iF l J i ' 2 (q,p-,p) + (q+p) a a B (q,p-,p) 


+ l -jY L ■ (q+p) A A{q,P',p) + ^y 11 • (q+p)A c (q,p*p) 


( 6 ) 


with 


BC (q,p-,p) = iL A (q,p;p)Yo+iZ c (q,p-,p)YL O) 


where q = q + u, p = p + u, and yll = f?y • u, y x = y — yll with 
m 2 = 1, 


T. F {q,p;p) 

A F (q,p;p) 


\[F(q 2 ,q4;p)+F{p 2 lP4 \p)\ , 
F(q 2 ,q 4 ',p) —F(p 2 ,p 4 ',p) 


with F = A,B,C of Eq. @. For comparison, we will also in¬ 
vestigate the so-called 1BC vertex, i.e., the first term of Eq. ([6]) 
alone. 

For the gluon propagator, we still employ the Landau gauge 
form with an infrared-dominant interaction modified by the 
chemical potential MM 


1 i k 2 - k2 +a>P 

g 2 (p)D pc (k,p) =4n~d-^e ^ 



kpk a \ 

k 2 J' 


( 8 ) 


The various parameters can be obtained by fitting meson prop¬ 
erties and chiral condensate in vacuum [391, 42] and we use 
CO = 0.5 GeV, d = 1 GeV 2 (with the rainbow approximation), 
d = 0.5 GeV 2 (with the BC vertex), d = 0.75 GeV 2 (with the 
1BC vertex), m u c i = 0, and m s = 0.115 GeV. The phenomeno¬ 
logical parameter a represents a reduction of the effective in¬ 
teraction with increasing chemical potential. This parameter 
cannot yet be fixed independently and its value will be dis¬ 
cussed in the following. 

The EOS of cold QM is obtained following Refs. ff4l[4Bl 
lir fl. All the relevant thermodynamical quantities at zero tem¬ 
perature can be computed from the quark propagator at fi¬ 
nite chemical potential, except a boundary value of the pres¬ 
sure, which is represented by a phenomenological bag con¬ 
stant Z?ds = 90 MeV fm . See Ref. lfi~4il for details. 


III. RESULTS AND DISCUSSION 
A. EOS of dense matter in beta equilibrium 

In order to study the structure of NS, we have to calculate 
the composition and the EOS of cold, neutrino-free, charge- 
neutral, and beta-stable matter characterized by two degrees 
of freedom pn and p e , the baryon and charge chemical poten¬ 
tials. The corresponding equations are 

Pi — biPB q iPe - ^ Piqi = 0 , 


(9) 
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FIG. 1. (Color online) Baryon density (upper panel) and pressure 
(lower panel) vs. the baryon chemical potential of NS matter for 
different models. The results for hadronic matter are shown with 
a black solid curve (N-BOB) and for quark matter with curves of 
different thickness for the rainbow approximation (RB-a), the 1BC 
vertex (lBC-a), and the BC vertex (BC-a), respectively. The verti¬ 
cal dotted lines indicate the positions of the phase transitions under 
the Maxwell construction. 


bj and qt denoting baryon number and charge of the particle 
species i = n,p,e,ji in the hadron phase and i = u,d,s,e,p 
in the quark phase, respectively. At low density we assume 
a hadron phase and will investigate a possible transition to 
quark matter at high density under the Maxwell or Gibbs con¬ 
struction. 

Fig. Q] shows the numerical results for hadronic matter and 
quark matter in beta equilibrium with different choices of 
the quark-gluon vertex and different parameters a given in 
Eq. ([8]). One can directly read off the phase transition between 
hadron matter and quark matter under the Maxwell construc¬ 
tion, which is presented as crossing point of the baryon and 
quark pressure curves (lower panel) and the projected verti¬ 
cal lines of the baryon density curves (upper panel). Though 
the DSM EOS is generally stiffer than that with the MIT bag 
model El, it strongly depends on the quark-gluon vertex and 
the parameter a. Within a given vertex ansatz, the pressure 
and density increase as a increases. For the same value of a, 
the pressure and density obtained with the 1BC or BC vertex 
are larger than those with the rainbow approximation. For the 
DSM with the rainbow approximation, we repeat the results 
with a = 1,2 10 , which were able to support hybrid stars 
with maximum masses larger than two solar masses. 



FIG. 2. (Color online) Effective bag constant of beta-stable quark 
matter as a function of baryon number density for different quark 
models. 


However, as stated above, a is a phenomenological param¬ 
eter that currently cannot be fixed independently. In the fol¬ 
lowing, we will not compare the results between different ver¬ 
tex ansatze with the same value of a, but instead we constrain 
the parameter a with different vertex ansatze by requiring the 
same phase transition point under the Maxwell construction. 
In other words, we employ a global reduction rate of the inter¬ 
action due to both the gluon propagator and the quark-gluon 
vertex by fitting the same phase transition point under the 
Maxwell construction. This requires that c<bc < G|isc < Gris- 
In particular, we find a = 0.8(1.7) for the 1BC ansatz and 
a = 0.25(0.93) for the BC ansatz corresponding to a = 1(2) 
for the rainbow approximation. With such a choice of a, we 
obtain results with the rainbow approximation and 1BC ansatz 
that are almost indistinguishable, while the pressure and den¬ 
sities with the BC ansatz are only slightly lower as functions 
of chemical potential. 

A useful quantity to characterize an effective quark model 
is the effective bag constant jus, 

fi(p B ) = £(p fl )-£f r ee(pB), (10) 

defined as the difference of the energy densities of beta-stable 
quark matter obtained for the DSM and free quark system, 
respectively. It is shown in Fig. [2] For all the three vertex 
ansatze, the effective bag constant is a monotonically increas¬ 
ing function of baryon density, presenting an important qual¬ 
itative difference from the generalized MIT bag model @], 
and are quantitatively much larger than the effective bag con¬ 
stants from the color dielectric model El and the Nambu- 
Jona-Lasinio model 1441. Rl. Once again, we find that also 
this quantity depends mainly on the rescaled a parameter, and 
not on the details of the vertex, in which case the rainbow ap¬ 
proximation and the 1BC ansatz are almost indistinguishable, 
while the results with the BC ansatz are only slightly higher 
at high densities. With larger reduction rate of the interaction 
the effective bag constant increase faster. 

We also investigate the phase transition under the more 
sophisticated Gibbs construction EUHSi, as detailed in 
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FIG. 3. (Color online) Pressure vs. baryon density of NS matter with 
the Gibbs (upper panel) and Maxwell (lower panel) phase transition 
construction for different quark models. 



FIG. 4. (Color online) Fractions of all particle species as functions of 
baryon number density with the Gibbs phase transition construction 
for different quark models. 


Eq. (17-22) of Ref. Q, which comprises a mixed phase re¬ 
gion of oppositely charged hadron and quark matter domains, 
while preserving the total charge neutrality. The Gibbs con¬ 
struction represents the zero-surface-tension limit of the cal¬ 
culations including finite-size effects (12i [0,0], whereas the 
Maxwell construction corresponds to very large surface ten¬ 
sion. The resulting EOSs p{ps) are shown in the upper panel 
of Fig. O while the results under the Maxwell construction are 
shown in the lower panel for comparison. Again one observes 
that the density range and the EOS of the mixed phase de¬ 
pend mainly on the global reduction rate, and are then almost 
indistinguishable with the three vertex ansatze. 

Fig. [4] shows more detailed information of the mixed phase, 
namely the fractions x- L = pi/P b of all species of particles as 
functions of baryon number density, where in the mixed phase 
with quark volume fraction % the average partial densities are 
given by p, = %p, for quarks and p, = (1 x)pi for hadrons. 
Also these quantities are almost the same for different vertex 
ansatze, once the properly scaled parameter a is chosen. 


B. Neutron star structure 


are all unstable and the corresponding results are not shown 
here. For the description of the NS crust, we have joined the 
hadronic EOS with the ones by Negele and Vautherin l5Qt] 
in the medium-density regime, and the ones by Fe ynm an- 
Metropolis-Teller lf5lll and Baym-Pethick-Sutherland |52| for 
the outer crust. 

The results are plotted in Fig. [5j where we display the grav¬ 
itational mass M as a function of the central baryon number 
density p c and the radius R. For all the three vertex ansatze 
and the chosen values of a, we obtain hybrid stars with the 
maximum mass lower than the pure nucleonic NS, but higher 
than two solar masses. From Figs. [4] and [5] one deduces that 
even in the most massive hybrid stars no pure quark phase ex¬ 
ists, but only a mixed phase of nuclear matter and QM in the 
inner core. Consistent with the results for the EOS, also the 
hybrid NS mass-radius relation depends mainly on the scaling 
parameter a and does then not allow to distinguish between 
the different choices of the quark-gluon vertex. To get a bet¬ 
ter understanding of the quark-gluon vertex in the DSM, one 
therefore needs refined and independent information about the 
parameter a, i.e., the in-medium modification of the gluon 
propagator. 


As usual, we assume that a NS is a spherically symmet¬ 
ric distribution of mass in hydrostatic equilibrium and obtain 
the stellar radius R and the gravitational mass M by the stan¬ 
dard process of solving the TOV equations i49il . We have 
used as input the EOSs with the Gibbs construction discussed 
above and shown in the upper panel of Fig. [3 The hybrid 
NS with phase transitions under the Maxwell construction 


IV. CONCLUSIONS 

We have further developed our DSM for QM under the BC 
and 1BC ansatze for the quark-gluon vertex in comparison 
with the previous rainbow approximation. We investigated the 
EOS of beta-stable quark matter in NS and the phase transi¬ 
tion between hadronic matter and QM in combination with a 
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FIG. 5. (Color online) Gravitational NS mass (M@ = 2 x 10 33 g) 
vs. radius (right panel) and normalized central baryon density (left 
panel; po = 0.17 fm~ 3 ) for different EOSs. 


microscopic hadronic EOS obtained within the BHF formal¬ 
ism. 

Due to the uncertainty of the in-medium gluon propagator 
in the DSM, which we currently approximate by a single pa¬ 


rameter a, we obtain a global modification of the interaction 
caused by both the quark-gluon vertex and the gluon propa¬ 
gator. The freedom of the parameter a allows then to fit the 
same phase transition point under the Maxwell construction 
for the different choices of the vertex. Consequently, we find 
that all observables, i.e., EOS, particle fractions, and stellar 
structure, mainly depend on the global reduction rate of the 
interaction, but are not sensitive to the particular ansatz for 
the quark-gluon vertex. With a proper choice of a and for 
all three vertex ansatze, the EOS from the DSM can support 
a two-solar-mass hybrid star, provided that the appearance of 
hyperons is excluded. 

These results focus the attention to the necessity of ob¬ 
taining model-independent information about the quark-gluon 
vertex and gluon propagator at finite baryon density. More in¬ 
vestigation on the quark-gluon vertex and gluon propagator 
from their own DSEs will be done in the future. 
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